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1. Introduction 

One distinct property of quantum entanglement from other classical correlations is 
its restricted sharability. For instance, if a pair of parties in a multipartite quantum 
system share maximal entanglement, then they can share neither entanglement [U [2] 
nor classical correlations [3 J with the rest. This is known as the Monogamy of 
Entanglement (MoE) J3], and it has been shown that this restricted sharability of 
quantum entanglement can be used as a resource to distribute a secret key which is 
secure against unauthorized parties [51 [6] . 

Whereas MoE is a restricted property of entanglement in multipartite quantum 
systems, the sharability itself is about the bipartite entanglements among the parties in 
multipartite systems. In other words, it is inevitable to have a proper way of quantifying 
bipartite entanglement for a good description of the monogamy nature in multipartite 
quantum systems. For this reason, certain criteria of bipartite entanglement measure 
were recently proposed for a good description of the monogamy nature of entanglement 
in multipartite quantum systems [TJ; that is, 

(i) Monotonicity: the property that ensures entanglement cannot be increased under 
local operations and classical communications. 

(ii) Separability: capability of distinguishing entanglement from separability. 

(iii) Monogamy: upper bound on a sum of bipartite entanglement measures thereby 
showing that bipartite sharing of entanglement is bounded. 

The first mathematical characterization of MoE was shown in three-qubit 
systems [1] using concurrence [8] as the bipartite entanglement measure. It is known 
as CKW -inequality named after its establishers, Coffman, Kundu and Wootters, and 
this CKW- type inequality was also shown for arbitrary multi-qubit systems later [2]. In 
other words, concurrence is a good entanglement measure for multi-qubit-systems that 
satisfies the criteria proposed in [7j. 

However, monogamy inequality using concurrence is know to fail in its 
generalization for higher-dimensional quantum systems [9j [7J. Furthermore, although 
MoE in multi-qubit systems is mathematically well-characterized in terms of 
concurrence, it is not generally true for other entanglement measures such as 
Entanglement of Formation (EoF) pi)]. In other words, MoE does not have CKW- 
type characterization in terms of EoF, and this exposes the importance of the choice 
of a bipartite entanglement measure to characterize MoE even in multi-qubit systems. 
Moreover, for possible generalization of monogamy inequality into higher-dimensional 
quantum systems, it is undoubtedly one of the most important and necessary tasks to 
have a proper way of quantifying bipartite entanglement. 

Renyi-a entropy pT| is a generalization of Shannon entropy [12] in terms of the 
non-negative real parameter a, and it has been widely used in the study of quantum 
information theory such as quantum entanglement and correlations [131 HH 03] • As 
a way to quantify the uncertainty of probability distribution, Renyi-a entropy shows 
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its selective characters with respect to a: The higher values of a, Renyi-a entropy is 
increasingly determined by consideration of only the highest probability events, whereas 
it increasingly weights all possible events more equally, regardless of their probabilities 
for lower values of a. For the case when a tends to 1, Renyi-a entropy converges to 
Shannon entropy 

Here, we show that the distinct character of Renyi-a entropy with respect to the 
parameter a also appears in establishing monogamy inequalities of entanglement in 
multipartite quantum systems. Although EoF that is based on Shannon entropy is 
known to fail for usual CKW-type characterization of MoE, Renyi-a entropy can still 
be shown to have CKW-type monogamy inequality for all case of a if it exceeds a certain 
threshold. We first provide an analytic formula for the bipartite entanglement measure 
based on Renyi-a entropy namely Renyi-a entanglement in two-qubit systems for a > 1, 
and we show that multi-qubit entanglement shows the usual CKW-type monogamy 
inequalities in terms of Renyi-a entanglement for a > 2. For a < 2, we conjecture with 
strong numerical evidence that Renyi-a entropy can provide a possible dual monogamy 
or polygamy inequality of multi-qubit entanglement for 0.83 — e < a < 1.43 + e with 
< e < 0.01 

This paper is organized as follows. In Section 12. 1[ we recall the definition of Renyi- 
a entropy, and Renyi-a entanglement for bipartite quantum states. In Section 12.21 we 
provide an analytic formula of Renyi-a entanglement for arbitrary two-qubit states. In 
Section 13.11 we derive a monogamy inequality of multi-qubit entanglement in terms of 
Renyi-a entanglement for a > 2, and we conjecture a polygamy inequality of multi-qubit 
entanglement using Renyi-a entropy with strong numerical evidences in Section 13.21 
Finally, we summarize our results in Section HI 

2. Renyi-a Entropy and Entanglement Measures 

2.1. Definition 

For a probability distribution P = {p^ where < Pi < 1 for all i and YliPi = 1> 
classical Renyi-a entropy is defined as 



for any positive a such that a^l. Throughout this paper, the logarithmic function is 
assumed to have base two otherwise specified. In the limiting case where a tends to 1, 
H a (P) converges to Shannon entropy, that is, 




(1) 




(2) 



For any quantum state p, its quantum Renyi-a entropy is defined as 



(3) 
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for any a > and a^l [TB]. For the quantum state p with its spectral decomposition 

p = Ei-MV'iXV'il, we have 

Sa(p) = H a (X), (4) 

where X = {Aj} is the spectrum of p, and thus, S a (p) converges to the von Neumann 
entropy of p for the case when a — > 1. In other words, Shannon entropy and von 
Neumann entropy are the singular points of classical and quantum Renyi entropies 
respectively, and those singularities are removable. For this reason, we will just consider 
Hi (P) = H (P) and Si(p) = S(p) for any probability distribution P and quantum state 
P- 

For a bipartite pure state \ip) AB , the von Neumann entropy of the reduced density 
matrix p A = ti b\^) ab{^\ ls known to be a good bipartite entanglement measure 

E{\^)ab) = S(p A ) = S(p B ). (5) 

With noticing that von Neumann entropy quantifies the uncertainty of the quantum 
state, this way of quantifying bipartite entanglement is based on the uncertainty of 
subsystem: More uncertainty on subsystems implies stronger quantum correlation 
between subsystems. 

A well-known way to generalize this concept of entanglement measure into mixed 
states is taking the minimum (or infimum) of the average entanglements 

E { {p AB ) =mm^2 Pi E(\tp) AB ) (6) 

i 

over all possible pure state decompositions of the mixed state pab = YliiVi^i) ab^^- 
This generalization is known as convex-roof extension, and Ei( P ab) is called the 
entanglement of formation of pab- 

As a generalization of EoF into the full spectrum of Renyi-a entropy [17], Renyi-a 
entanglement of a bipartite pure state |V')^4.b is defined as 

E a m AB ) = S a ( PA ), (7) 

where pa = tr 5 |^) AB (^|, and for a mixed state Pab> its Renyi-a entanglement is defined 

as, 

E a {pab) = min^2 PiE a (\ipi) AB ), (8) 

i 

where the minimum is taken over all possible pure state decompositions of pab — 
YliViH'i) AB^il- Similar to EoF for bipartite quantum states, Renyi-a entanglement 
is defined based on the uncertainty of subsystems, which has EoF as a special case when 
a — > 1. 

It is direct to check that E a (pab) = if and only if pab is a separable state, and 
furthermore, Renyi-a entanglement is also known as entanglement monotone: it is not 
increased under local operations and classical communications. 
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2.2. Analytical formula for two-qubit systems 

Let us recall the definition of concurrence. For any bipartite pure state I^/abi 
concurrence, C(\(f>) AB ) is defined as [8] 



C(\4>)ab)= y/2(l-tTp> A ), (9) 
where pa = tr B (|0) AB {4>\), and for any mixed state Pab, its concurrence is defined as 
C(p A b) = min^2p k C(\(f) k } AB ), (10) 

k 

where the minimum is taken over all possible pure state decompositions, pab = 

For any two-qubit mixed state pab in B (C 2 <8> C 2 ) , its concurrence is known to 
have an analytic formula [8], that is, 

C(pab) = max{0, Ai - A 2 - A 3 - A 4 }, (11) 

where Aj's are the eigenvalues, in decreasing order, of v/ y 'Pab 'Pab v 'Pab and Pab = 
a y ® o- y p* AB Oy ® a y with the Pauli operator a y . Furthermore, the relation between 
concurrence and EoF of a two-qubit mixed state pab (or a pure state \4>) A b e C 2 (S>C d ), 
can be given as a monotone increasing, convex function 8 [8], such that 

Ef(pAB) = £(Cab), (12) 

where 

£{x) = H(- + -Vl-x 2 ^), for0<ac<l, (13) 

with the binary entropy function H(t) = —tlogt — (1 — t) log(l — t). In other words, 
the analytic formula of concurrence as well as its functional relation with EoF lead us 
to an analytic formula of EoF for two-qubit states. 

For any two-qubit pure state (or any pure state with Schmidt-rank less than or 
equal to two) with its Schmidt decomposition |V ; )ab = V^o|00)ab + V^i|H)ab> ^ s 
Renyi-a entanglement is 

E a{\lp) A B) = S *(PA) 

= -±_log(\« + \«). (14) 
Now, for each a > 0, by defining an analytic function 



fa{x) ■= T-^— log 

1 — a 



X 



(15) 



on < x < 1, it can be directly checked that 

Ea(mAB)=fa(C(\lP) AB )), (16) 

where C(\i/j) ab ) is the concurrence of \ip) AB - Thus, for each a > 0, we have a functional 
relation between the concurrence and Renyi-a entanglement for two-qubit pure state 
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\4>) AB - Here, we note that for the case when a tends to 1, f a (x) converges to the 
function £(x) in Eq. (ITBl ; that is, 

lim f a (x)=S(x). (17) 

a— >1 

For a mixed state pab, we have the following theorem. 
Theorem 1. For each a > 0, if 



f a {x) = — ^— log 
1 — a 



x 



is a monotonically increasing and convex function, then 

E a (PAB) = fa {C{PAB)) , (19) 

for any two-qubit state pab where C(pab) is the concurrence and E a (pab) is the Renyi-a 
entanglement of pab ■ 

Proof. Suppose pab = YliiVi^i) ab(^I is an optimal decomposition for E a ( P ab) such 
that E a (pab) = J2iPi E * (I^)ab)' then 

Ea(pAB) = ^2piE a (\^i) AB ) 

i 

= J2Pifa(C(\lpi) AB )) 

i 

AB) 

>fa(C(p AB )) (20) 

where the second equality is by the relation between concurrence and Renyi-a 
entanglement for pure states in Eq. ffl6l) . the first inequality is by the convexity of 
f a , and the last inequality is by the monotonicity of f a and the definition of C (pab)- 

Due to the analytic formula of concurrence for two-qubit states [8], we can always 
assume the existence of an optimal decomposition pab = Y^iPilrfi) AB^il sucn ^ na ^ 

C(pab) = Y,P^ C ^)ab) (21) 

i 

and 

C(m AB )=C( PAB ) (22) 
for all i. Now, we have 

f a {C (pab)) = fa {^2 PiC (I^)ab) 

= ^PifaiC(\^i) AB )) 

i 

= ^PiEad^i)^) 
i 

> E a {p AB ) (23) 
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where the second equality is by Eq. 
E a (pab)- 

Thus, by Eqs. (1201) and ( |23l) . we have 

E a (pab) = fa {C{pab)) , 
which completes the proof. 



and the inequality is by the definition of 

(24) 
□ 



Theorem [T] implies that for any positive a such that f a {x) in Eq. (fT5|) is 
monotonically increasing and convex, the analytic formula of concurrence for a two- 
qubit state pab in Eq. ffTTj) together with the functional relation between C (pab) and 
E a {.Pab) in Eq. (fT9l) provide us an analytic formula for Renyi-a entanglement of pab- 
The rest of this section is mainly contributed to the analytic proof for the range of a 
where f a is monotonically increasing and convex. 



Theorem 2. For any real a > 1, 



fa (x) 



1 



a 



log 



1-vT 



i + vT 



(25) 



is a monotonically increasing and convex function for < x < 1 . Furthermore, the 
monotonicity and convexity of f a (x) are strict in the sense that f a (x) < f a (x') and 
f a (\x + (1 — \)y) < \f a (x) + (1 — X)fa(y) for any < A < 1 and < x, x', y < 1 such 
that x < x' . 

Proof. For the case when a tends to 1, f a {x) converges to S(x) in Eq. (fl~3l) which is 
already known to be monotonically increasing and convex [8]. Here, we consider the 
case where a > 1. 

Let us define a function 



9a[X) 



log 



X/ 



then we have 



fa{x) 



1 



-9a{x) + 



a 



(26) 



(27) 



a — 1 a — 1 

which implies that the convexity and the monotonicity of f a {x) follow from those of 
g a {x) for a > 1. Furthermore, g a {x) is an analytic function on < x < 1, thus the 
convexity and monotonicity of (^(x) can be assured by showing its first and second 
derivatives are nonnegativity. 

By taking the first derivative of g a {x), we have 



dg a (x) 
dx 



ax 



(l + VT^Y 1 - (l - VT^) 



a-l 



(l - VT^) + (l + vT 



> 



for all a > 1 and the equality holds only at the boundary, that is x = or x 
other words, g a {x) is a strictly monotone-increasing function for < x < 1. 



(28) 



1. In 
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For the second derivative of g a (x), we have 

d 2 9a(x) _ B {At 1 - AT 1 ) {A\ + A%) 



dx 2 



+ B 



2 (AT 1 - Af 1 ) 2 - 4(a - 1) 



x 



„2a-2 



X 



(29) 



where A x = 1 + y/l - x 2 , A 2 = 1 - v 7 ! - £ 2 and 5 = a/ [(1 - x 2 ) (Af + A%f] . The 
binomial series for A"" 1 and Af 1 lead us to 

ct-1 



Af 1 = 1 + vT 



l + (a-l)Vl-x 2 + 



and 



^r 1 



a-1 



1 - [a- l)Vl -x 2 + 



(g-l)(g-2) 
2! 



(a-l)(ot-2) 
2! 



(30) 



X z - • 



(31) 



and thus 



Af 1 _ Af 1 = 2(a - l)VT^r2 + Ci > 2 ( a _ 1)VT^ 



x z 



1 + Af 1 = 2 + C 2 > 2, 

for some non-negative C\ and C^- Now, let 

(Af 1 - Af x ) (A° + A2) 



r 2 = x 2 (at 1 -at 1 ) 2 



vT 



if- 



then by Eq. (|32|) . we have 

r a > 4(a - 1), T 2 > 4(a - 1) V(l - x 2 

and 



(32) 



(33) 



(34) 



d;2 J B[r 1 + r 2 -4(«-i)x 2a - 2 ] 

> B [4(a - 1) - 4(a - 1)x 2q " 2 + 4(a - 1) V(l - x 2 

> 0, 



(35) 



for any a > 1 where the equality holds if and only if x = 1. Thus, the first and second 
derivatives of g a (x) are nonnegative for < x < 1 and strictly positive for < x < 1, 
which implies the strict monotonicity and convexity of f a (x) on < x < 1 for a > 1 . □ 

Now, we have the following corollary, which is the first main result of this paper. 

Corollary 1. For any a > 1, a two-qubit state Pab has an analytic formula for 
its Renyi-a entanglement such that E a (pab) = fa (C(pab)) where C(pab) is the 
concurrence of pab, and 



fa(x) 



1 — a 



log 



+ 



(36) 
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x 



(a) (b) (C) 

Figure 1. d for a — 0.83 are illustrated in (a) and (b) on the domain < x < 1 
and 0.99 < x < 1 respectively. The convexity of f a (x) is violated for a = 0.82 in (c) 
where d has negative values for x around 1. In Picture (b), the vertical axis is 

scaled by 10~ 2 . 



In fact, %M i n Eq. (EBT) can be easily checked to be nonpositive for < a < 1, 



Ax 

and together with Eq. (TJTD, we have ^™ = ^-^^ > for < a < 1. In other 
words, the function f a (x) in Eq. ffTBT) is a monotone-increasing function on the domain 
of < x < 1 for any positive a. 

For the convexity of f a (x) with < a < 1, we first note that the continuity of 
d j^i^ in Eq. (129]) with respect to a assures the positivity of d for a slightly less 
than 1. However, for general a between and 1, we tried a numerical way of calculations 
for the second derivative of f a (x) which is illustrated in Figured! 

We have tested d g ^f^ for various values of a between and 1 and observed that 
d is nonnegative for a > 0.83 (Figure [1] (a) and (b)), whereas its positivity 
is violated for a = 0.82 with x around 1 (Figure [T] (c)). In other words, f a (x) is 
numerically observed to be still convex for 0.83 — e < a < 1 with some positive e such 
that < e < 0.01. Thus, together with Theorem [2j we have the following conjecture. 



Conjecture 1. For any positive a > 0.83— e with some positive e such that < e < 0.01, 

(37) 



fa{x) = — ^— log 

1 — a 

is a convex function for < x < 1. Furthermore, for this range of a, any two-qubit 
state pab has an analytic formula for its Renyi-a entanglement such that E a (pab) = 
fa (C(pab)) where C(pab) is the concurrence of pab- 

3. Entanglement Constraint in Multi-party Quantum Systems 



In this section, we establish a mathematical formulation for the monogamous and 
polygamous properties of multi-qubit entanglement in terms of Renyi-a entanglement. 
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First, we show an important property of the function f a (x) in Eq. (1151) . By using the 
property of f a (x) as well as the functional relation between Renyi-a entanglement and 
concurrence for two-qubit states in previous section, we derive a monogamy inequality 
of multi-qubit entanglement in terms of Renyi-a entanglement for a > 2. We also 
conjecture a polygamy inequality of multi-qubit entanglement in terms of Renyi-a 
entanglement for a around 1. 

3.1. Monogamy of multi-qubit entanglement 

For a three-qubit pure state \tp) ABC CKW-inequality was shown as [I], 

Ca(bc) — Cab + ^ac> (38) 

where Ca(bc) — Cd^P) a(bc)) * s ^ ne concurrence of a 3-qubit state |'0)a(bc) w ^h respect 
to the bipartite cut between A and BC, and Cab and Cac are the concurrences of the 
reduced density matrices onto subsystems AB and AC respectively. Later, Eq. (I38l was 
generalized into arbitrary multi-qubit systems as 

^A^Aa-An) ^ C MA 2 H H ^-AiA n 5 (39) 

for an n-qubit state /?Ai-A n 0- 

In this section, we provide a class of monogamy inequalities in multi-qubit systems 
in terms of Renyi-a entanglement. Before we prove this, we first show a property of the 
function f a (x) that plays a crucial role in the proof of Renyi-a entanglement monogamy. 

Theorem 3. For any real a > 2 and the function, 



f a (x) = log 

1 — a 



(40) 



we have 

f a (Vx^+f) > f a {x) + f a {y) (41) 
for < x, y < 1 such that < x 2 + y 2 < 1. 

Proof. Let P = {(x, y)| < x, y, x 2 + < 1} and h a (x,y) be a function defined on 
the domain D such that 

h a (x,y) := f a (V* 2 + y 2 ) ~ fa(x) - fa{y)- (42) 

Then, it is enough to show that h a (x,y) is a non-negative function on D. 

Because h a (x, y) is continuous on the domain T> and analytic in the interior of T>, its 
maximum or minimum arises at the critical points or boundary of T>. Here, we will show 
that h a (x,y) has neither vanishing gradient in the interior of T> nor negative function 
values on the boundary of T>, which implies that h Q (x,y) is a non-negative function on 
the domain T>. 

First, let us consider the gradient of h a (x, y), 

V Wl , j)= (^),^!) (43) 
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dh a (x,y) 
dx 



Cx 


(i + v 


l\ — x 2 — y 2 


a-l 

-( 


1 — \/l — x 2 — y 2 


a-l 




a/1 - x 2 - y 2 


(l-^l-x 2 -y 2 } 


+ (l + ^/l-x 2 


-y 2 


)"] 



Cx 



(l + VT^) a 1 - (l - VT^x 2 ) 



a-l 



v 7 !^ 1 (i - v 7 ! 3 ^) + (i + VT^y 



(44) 



and 



dh a (x,y) 
dy 



Cy 


(i + v 


/l — x 2 — y 2 


a-l 


1 — a/1 — x 2 — y 2 


a-l 




a/1 — x 2 — y 2 


(l-^l-x 2 - y 2 J 


+ (l + ^1-x 2 


-y 2 





Cy 



'i-r 



a-l 



i- Vi -r 



a-l 



'1-y 2 



1 



l-y 2 + [l + Jl-y 



(45) 



with C = Now, suppose V/i a (xi, yi) = (0,0) for some (#1,2/1) in the interior of T>, 
that is, < xi, yi < 1 and < x 2 + y 2 < 1. Because both xi and are nonzero, from 
Eqs. (J44J) and we have 

(1 + v^^)- 1 _ (i _ vr^?)^ 1 " 



vr^? [(i - a/t^T + (1 + v/r^?)"] 







a-l 

- 


(l - Vl-2/i 2 ; 


a-l 




-2/i 2 


[(i - 


-yi 2 ) 


+ (i + 


■ yi 2 ^ 


n 



(46) 



Furthermore, by defining a function l a (x) on < x < 1 such that 

, «— 1 /_ /- a—1 



l a (x) 



(i + v / r^) - (i - v / r^ 1 ) 



v 7 !^ 2 (l - VT^) + (l + VT^ 2 ) 



Eq. (T46|) can be rewritten as 
l a (xi) = l a {yi). 
Here, we note that the first derivative of l n (x) is 



(47) 



(4f 



dZ tt (g) 
dx 



5 



(A? + A«) (Ar 1 - At 1 ) 



+ (Ar 1 -A%- i y 



x/r^ 2 



4S (a - 1) x 



2a- 4 



(49) 
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with A\ — 1 + vl — x 2 , Ao = 1 — \/l — x 2 and 5 = and from the 

(l- x a)(Af+^) ! " 



> 4(a - 1) 



inequalities in Eq. (1321) . we have 

(A? + A%) (AT 1 -AT 1 ) 



X 



2 



[AT 1 -AT 1 ) 2 >4(a-l) 2 (l-x 2 ), (50) 

which implies 

^4^- > 45 [(<* - 1) + (" - 1) 2 (1 - x 2 ) - (a - 1) x 2Q - 4 l . (51) 
ax L J 

For the region of a > 2 and 0<x<l,a — lin the second term of the right-hand 

side of Eq. (!5T!) is always larger than or equal to 1, whereas x 2ct ~ 4 in the last term is 

always less than or equal to 1, thus 

^4^- > 45 \(a -l) + {a- 1) 2 (1 - x 2 ) - (a - 1) x 2a ~ 4 ] 
ax L J 

> AB [(a - 1) + (a - 1)(1 - a; 2 ) - (a - 1)] 
= 4B(a-l)(l-x 2 ) 

> 0, (52) 

and the last inequality is strict for < x < 1. In other words, / a (x) is a strictly 
monotone- increasing function on < x < 1, and therefore Eq. (T4*8]) implies Xi = yi. 
(Also note that the possible range of x\ now becomes < x\ < -^=.) However, from 
Eqs. f l4"4"|) and ([15]) . V/i a (xi, yi) = (0,0) together with Xi = |/i implies that 

=Z a (a; 1 ), (53) 

which contradicts to the strict monotonicity of l a (x) for < x < 1. Thus, V/i Q (xi, y±) ^ 
(0, 0) for any (xi, in the interior of T>. 

Now, let us consider the function value of h a (x,y) on the boundary of T>, that is 
x = or y = 0, or x 2 + y 2 = 1. If one of x or y is 0, then it is clear to check that 
h a (x, y) = 0. For the case when x 2 + y 2 = 1, we have 

h a (x,y) = f a (yx 2 + y 2 ^j - f a (x) - f a (y) 

= fa(l)-fa(x)-f a (VT^) 

2a - + ^-log[(l + x 2 ) Q +(l-x 2 ) Q ] 
1 a — 1 

(i - v / r 3 ^ 2 ) a + (i + v 7 ! 3 ^ 2 )" 

=: m a (x) (54) 

where m a (x) is a one-parameter real- valued analytic function defined on < x < 1. For 
a > 2, it is also straightforward to check that the function m a (x) has only one critical 
point at x — ^= through the domain < x < 1. Furthermore, all the function values of 
m a (x) at the critical point and the boundary are nonnegative, that is, we have 

h a (x,y) = m a (x) > (55) 



a 
1 . 

+ 7 log 

a — 1 



Monogamy and polygamy for multi-qubit entanglement using Renyi entropy 13 

for x 2 + y 2 = 1 and < x < 1. Thus, h a (x, y) is a non- negative function on the domain 
"D, and this completes the proof. □ 

Now, by using Theorem [3] together with Corollary [1] in previous section, we have 
the following theorem, which is the secondary result of this paper. 

Theorem 4. For a > 2 and any multi-qubit state PA^A 2 -A n ^ n & ^C 20 "^, we have 

E a (pAi(A 2 -A n )) > E a (pAiAa) H h E a (p Al A n ) , (56) 

where E a (pAi(A 2 -A n )) is the Renyi-a entanglement of p Al A 2 -A n with respect to the 
bipartite cut between A\ and the others, and E a (p AlAt ) is the Renyi-a entanglement 
of the reduced density matrix pA 1 A i on two-qubit subsystem A\Ai for i = 2, . . . , n. 

Proof. From Eq. fl39l) . we have 



C Al (A 2 -A n ) > \]C Al A 2 H ^ C MA^ ( 57 ) 

where Ca 1 (a 2 -a„) and Cam* are the concurrences of p Al ( A2 - An ) and pA 1 A i respectively. 
Now, let Pa^-Au) = J2iPi\4>i) Ai(A 2 -A n )(V , i| be an optimal decomposition for 

E a (p Al (A 2 -A n )) such that E a (p Al (A 2 -A n )) = J2iPi E a (\ipi} Al (A 2 -A n )) ■ Because each 
A 1 (A 2 -A n ) ^ n ^ ne decomposition has Schmidt-rank less than or equal to two (they are 
2 ® d pure states for d = 2® n_1 ), its concurrence and Renyi-a entanglement are related 
by the function f a (x) in Eq. ffT9l . that is, 

E « (\^)A l{ A 2 ...A n) ) = U (cm M A 2 ...A n) )) (58) 
for each i. Thus, we have 



E a {PA 1 (A 2 -A n )) = ^PiE a (\ll>i) Al (A 2 -A n ) 
i 

= YjPrf* ( c (i^)a i( a 2 ...a„) 



> fa \^p£(\A) Al{M ... An) ) 

> fa (C Al (A 2 -A n )) , (59) 

where the first inequality is by the convexity of f a (x), and the second inequality is by 
the definition of concurrence and the monotonicity of f a (x). Furthermore, by Eq. ([5j 
together with Theorem [3], we have 

fa {C Al ( A2 - An )) > fa (y\JC- 2 Al A 2 "I f" C A lAn 

> fa (C Al A 2 ) + fa [\l C Al A A H ^ C Al A n ) 



> fa {CaiA 2 ) H 1" fa {CaiA„) 

= E a (p Al A 2 ) H \~ E a (p AlAn ) , 



(60) 
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where the first inequality is by the monotonicity of f a (x), the other inequalities are by 
iterative use of Theorem [31 and the last equality is by the functional relation of Renyi-a 
entanglement and concurrence for two-qubit states. Thus, Eqs. ( 1591) together with ( |60j) 
complete the proof. □ 

In fact, it was recently shown that Renyi-a entanglement for a = 2 can be 
used to establish a monogamy inequality for multi-qubit systems by straightforward 
calculation [18]. However, Theorem H] says that the monogamous property of 
entanglement in multi-qubit systems can be mathematically characterized in terms of 
Renyi-a entanglement for all positive real number a larger than or equal to 2. 



3.2. Polygamy of multi-qubit entanglement 

In previous section, we have established the monogamy inequalities of multi-qubit 
entanglement in terms of Renyi-a entanglement for all positive real number a > 2. 
Here, we consider the case when < a < 2, and claim another kind of entanglement 
constraint in multi-qubit systems using Renyi-a entropy. 

Let us first recall the definition of Entanglement of Assistance (EoA) [19] for a 
bipartite state pab, that is, 

E a ( PAB ) = maxJ2PkE(M AB ) , (61) 

k 

where the maximum is taken over all possible pure state decompositions of pab — 
J2kPk\i J k) Asi^kl- Here, we note that EoA in Eq. ( l6Ti) is clearly a mathematical dual 
to EoF in Eq. because one of them is the maximum average entanglement over 
all possible pure state decompositions whereas the other is the minimum. Moreover, 
by introducing a third party C that has the purification of pab, E a (pAs) can also 
be considered as the maximum achievable entanglement between A and B assisted by 
C. (This is the reason why it is called the assistance.) In other words, E a (pAB) is 
the maximal entanglement that can be distributed between A and B assisted by the 
environment; therefore, EoA is also physically dual to the concept of formation. 

Similar to the duality between EoF and EoA, we can also have a dual concept to 
concurrence: Concurrence of Assistance (CoA) [20] for a bipartite state pab is defined 
as 

C a {pAB) = max^2p k C(\^ k ) AB ), (62) 

k 

where the maximum is taken over all possible pure state decompositions of pab = 
YlkP^k) AB(ipk\- Furthermore, it was shown that there exists a different kind of 
entanglement constraint in multi-qubit systems in terms of CoA [21]. More precisely, 
for any pure state \i>) a^-a in an n-qubit system, we have 

Cl { A 2 ...A n) <(C a A 1 A 2 ) 2 + --- + (C a Al Aj 2 , (63) 
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where C Al ( A2 ... An ) * s the concurrence of \i>) Al .., A with respect to the bipartite cut 
between A 1 and A 2 ■ ■ ■ A n , and C AlA . is the CoA of the reduced density matrix p Al Ai for 
i = 2, . . . , n. 

Eq. (163|) is known as the dual monogamy or polygamy inequality of entanglement in 
multi-qubit systems: Whereas concurrence can be used to characterize the monogamy 
of multipartite entanglement as in Eq. fl39l) . its dual quantity, CoA can also be used for 
the dual monogamy of multipartite entanglement. Later, it was shown that polygamy 
of multi-qubit entanglement can also be characterized in terms of EoA, that is, for any 
multi-qubit pure state A ... A , 



E (J^ i( a 2 ...a„)J < E a ( PAl A 2 ) + ■■■ + E a (p AlAn ), (64) 

where E Al ( A2 ... An )) — S (p A ) is the entanglement of \ip) Al ... An with respect to the 
bipartite cut A\ and A 2 ■ ■ ■ A n , and E a (p AlAi ) is the EoA of the reduced density matrix 
PAxAi for i — 2, • - • , n [22]. More recently, a general polygamy inequality of multipartite 
entanglement in arbitrary dimensional quantum systems was proposed by using an 
analytical upper bound of CoA [23] . 

Here, we claim that the polygamous property of multi-qubit entanglement can also 
be shown by using Renyi-a entropy for a around 1. In fact, we can intuitively expect 
this: Due to the continuity of Renyi-a entropy with respect to a, we can easily expect 
that the inequality in Eq. ( |64j) would also be true by using Renyi-a entropy, instead of 
von Neumann entropy, where a is in a region of small perturbation from 1. However, we 
conjecture, with strong numerical evidences, a specific region of a where the polygamy 
inequality can be obtained using Renyi-a entropy. 

For any bipartite state pab, let us first define a dual quantity to Renyi-a 
entanglement as 

K(Pab) ■= m&x^2p k E a (\ip k ) AB ) , ( 65 ) 
k 

where the maximum is taken over all possible pure state decompositions of p A s = 
YlkPkl^k) AB^fcl' ana - can it R^nyi-a Entanglement of Assistance (REoA). 

Similar to the functional relation between concurrence and Renyi-a entanglement 
for two-qubit states in Eq. (fl9|) . the same function f a (x) can also relate REoA of a 
two-qubit state with its CoA. For a two-qubit state pab, let pab = X^.PilV'i) ab^I De 
an optimal decomposition for its CoA, that is, 

C a (p^) = $>C(|^U), (66) 

i 

where C a (pab) is the CoA of pab defined in Eq. ( 1621) . By Theorem [2], f a (x) is convex 
for the range of a > 1, thus, for this range of a, we have 



'ABJ 



f a (c a ( PAB )) =f a P[>c(hM. 

< Ys^^^Ab)) 
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(a) Graph of hi 9 (x,y) for 
x 2 + y 2 < 0.008. 




(b) Graph of hiAi(x,y). 




1.5 



(c) Graph of hi.&{x,y). 



Figure 2. (a), (b) and (c) illustrate the graphs of h a (x,y) for a = 1.9, 1.44 and 1.43 
respectively. In each picture, the variables x and y are reparameterized as x = C cost 
and y = Csint for < C < 1 and < t < tt/2. The vertical axes in (a), (b) and (c) 
are scaled by 10~ 12 , 10~ 3 and 10 -3 respectively. 

= ^2piE a {\i)i) AB ) 

i 

< K(Pab), (67) 

where the first inequality is by the convexity of f a (x) and the last inequality is by the 
definition of REoA. Furthermore, based on Conjecture [U we also claim that Eq. fIBTl) is 
true for a > 0.83 — e with small e. 

Now, let us consider a property of h a (x,y) in Eq. f H2i) . From Theorem [3], we first 
note that the continuity of h a (x, y) in Eq. (142]) with respect to a assures its positivity for 
a slightly less than 2. In other words, the monogamy inequality in Eq. (1561) of TheoremH] 
that is analytically proven for a > 2 is also true for this region of a. However, h a (x, y) 
has negative function values for most case of < a < 2, which is illustrated in Figure [2j 
In Figure [2] (a), negative function values of h a (x,y) are observed near (x,y) = (0,0) 
when a = 1.9. Furthermore, the region in the domain of h a (x, y) with negative function 
values is getting larger as a decreases. (Figure [2] (b)) 

In fact, for the limiting case of a — > 1, h a (x,y) was analytically shown to be a 
non-positive function [22], that is, 

£(Vx^+^)<£(x)+£(y), (68) 

for non-negative x and y such that < x 2 + y 2 < 1 where 8{x) = liniQ,^! f a (x) is defined 
in Eq. ( TT31) . Again, due to the continuity of h a (x,y) with respect to a, we can assure 
that h a (x,y) is nonpositive everywhere when a is around 1. For a specific region of 
a, we have numerically tested for various values of a, and h a (x,y) is observed to have 
non-positive function values for < a < 1.43 + e with < e < 0.01, which is illustrated 
in Figure [2] (c). 
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Conjecture 2. For any positive < a < 1.43 + e with some positive e such that 
< e < 0.01 and the function 



fa( X ) = z log 

1 — a 



(69) 



defined on the domain V = {(x, y)\0 < x, y, x 2 + y 2 < 1}, we have 

fa(V^+¥) < fa(x) + f a (y). (70) 

Now, Conjectures [Hand [2] lead us to the following theorem, which is the last main 
result of this paper. 

Theorem 5. For a multi-qubit pure state \ip) A ...a an< ^ an y rea ^ a suc ^ that 0.83 — e < 
a < 1.43 + e with < e < 0.01, we have 

Ea (\i>) Al (A 2 -A n) ) < K(PA lA2 ) + ■■■ + E a a ( PAlAn ), (71) 

where E a a 1 (a 2 --a„)) = $a (pa) is the Renyi-a entanglement of \4>) a 1 (a 2 ---a u ) 
respect to the bipartite cut A\ and A 2 ■ ■ • A n , and E®(p AlA .) is the REoA of the reduced 
density matrix PA 1 A i for i = 2, • • • , n. 

Proof. The proof method follows the construction used in |22j. From the polygamy 
inequality of entanglement in multi-qubit systems in Eq. (163]) gives us an inequality 



e Ali A 2 -A n) < y^W 2 + • • • + (c% lAn y. (72) 

First, let us assume that {C AlA2 ) 2 + • • • + {C\ iAn ) 2 < 1, then we have 

E <* (|^)a 1 (A 2 ... j 4„)J = M C Ai(A 2 -A n )) 



<f a (\/(C a A l A 2 ) 2 + --- + (C a A 1 



\n) 2 



< fa {C a AlA2 ) + fa (^(CW 2 + " " " + 
— fa ipAtA^) + fa (Q1A3) + ' ' ' + fa if-AxAn) 

<E a a (p Al A 2 ) + --- + E a a (p Al A n ), (73) 

where the first equality is by the functional relation between the concurrence and the 
Renyi-a entanglement for 2®d pure states, the first inequality is due to the monotonicity 
of the function f a (x), the second, third and forth inequalities are obtained by iterative 
use of Eq. flTOj) in Conjecture [2], and the last inequality is by Eq. fl67|) . 

Now, assume that (C a AlM ) 2 + ••■ + (C% lA J 2 > 1. Since E a (\*P) Al (A 2 -A n) ) = 
'S'a(PAi) < 1 for any multi-qubit pure state \ip) a 1 (a 2 -a„)i ^ * s enou gh t° show that 
Ea(.PA 1 A 2 ) + • — I- E^,(pa 1 a„) > 1- Here, we note that there exist k e {2, . . . , n — 1} that 
satisfies 

(a*) 2 + ■ • • + ( C W 2 < 1. ( Ca MA 2 ) 2 + ■■■ + (C a AlAk+ f > I- (74) 
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By letting 

T := {C a MM f + ■ ■ ■ + (C a AlAk+1 f ~ 1, (75) 

we have 

1 = fa U(C a Al A 2 ) 2 + --- + (Cl lAk+1 ) 2 -T) 

< u (/(aj 2 +-+(a4) 2 ) + /« (\A C W 2 - T ) 

<^(p AlA2 ) + --- + ^(p AM J, (76) 

where the first inequality is by using Eq. ( 1701) with respect to (C AiA2 ) 2 + • ■ ■ + (C AlAk ) 2 
and (C AlA ) 2 — T, the second inequality is by iterative use of Eq. (jTDjl on (C^ iA J 2 + 
h {C% a ) 2 , and the last inequality is by Eq. ( ISTj) . □ 

4. Conclusion 

By using Renyi-a entropy, we have established a class of monogamy and polygamy 
inequalities of multi-qubit entanglement. We have shown that monogamy of multi-qubit 
entanglement can have CKW-type characterization in terms of Renyi-a entanglement 
for a > 2, and conjectured the possible region of a for polygamy inequality in terms 
of REoA. Although the specific region of a for polygamy inequality is supported by 
numerical evidences, the existence of such region is based on the continuity of Renyi-a 
entropy, which is analytically provable. Thus, the conjecture we claimed here is highly 
reasonable because it has both numerical and analytical reasons. 

Multipartite entanglement is known to have many inequivalent classes, which 
are not convertible to each other under Stochastic Local operations and classical 
communications (SLOCC) [21]. Furthermore, the number of inequivalent classes 
increases dramatically as the number of parties increase [25]. Not like bipartite 
entanglement, the existence of inequivalent classes of multipartite entanglement implies 
that the states from different classes are hardly comparable to each other in such a way 
of comparing a single parameter that quantifies their entanglement. This is one of the 
main difficulties in the study of multipartite entanglement. 

Whereas the interconvertibility of quantum states under SLOCC gives us an 
operational way to classify multipartite entanglement, entangled states from different 
classes can also reveal different characters with respect to their monogamy and polygamy 
properties. For example, three-qubit systems are known to have two inequivalent classes 
of genuine three-qubit entanglement, the Greenberger-Horne-Zeilinger (GHZ) class [26J 
and the W-class |24j . In terms of monogamy and polygamy relations, CKW and its 
dual inequalities are saturated by W-class states, while the differences between terms 
in the inequalities can assume their largest values for GHZ-class states. In other words, 
monogamy and polygamy of multipartite entanglement can also be used for an analytical 
characterization of entanglement in multipartite quantum systems. 
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The class of monogamy and polygamy inequalities of multi-qubit entanglement we 
provided here consists of infinitely many inequalities parameterized by a, and each 
inequality is based on the distinct character of Renyi-a entropy with respect to a. We 
believe that this selective choice of our monogamy and polygamy inequalities will leads 
us to an efficient way of analytic classification of multi-qubit entanglement. Moreover, 
our result will also provide useful tools and strong candidates for general monogamy 
and polygamy relations of entanglement in multipartite higher-dimensional quantum 
systems, which is one of the most important and necessary topics in the study of 
multipartite quantum entanglement. 
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